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In this series of lectures I discuss the basic principles and the modelling of the chemical evolu- 
tion of galaxies. In particular, I present models for the chemical evolution of the Milky Way 
galaxy and compare them with the available observational data. From this comparison one can 
infer important constraints on the mechanism of formation of the Milky Way as well as on stel- 
lar nucleosynthesis and supernova progenitors. Models for the chemical evolution of elliptical 
galaxies are also shown in the framework of the two competing scenarios for galaxy formation: 
monolithic and hierachical. The evolution of dwarf starbursting galaxies is also presented and 
the connection of these objects with Damped Lyman- a systems is briefly discussed. The roles 
of supernovae of different type (I, II) is discussed in the general framework of galactic evolu- 
tion and in connection with the interpretation of high redshift objects. Finally, the chemical 
enrichment of the intracluster medium as due mainly to ellipticals and SO galaxies is discussed. 



1. Basic parameters of chemical evolution 

Galactic chemical evolution is the study of the evolution in time and space of the 
abundances of the chemical elements in the interstellar gas in galaxies. This process is 
influenced by many parameters such as the initial conditions, the star formation and 
evolution, the nucleosynthesis and possible gas flows. 

Here I describe each one separately: 

• Initial conditions- One can assume that all the initial gas out of which the galaxy 
will form is already present when the star formation process starts or that the gas is 
slowly accumulated in time. Then one can assume that the initial chemical composition 
of this gas is primordial (no metals) or that some pre- enrichment has already taken 
place (e.g. Population III stars). As we will see in the following, different assumptions 
are required for different galaxies. 

• The birthrate function- Stars form and die continuously in galaxies, therefore a 
recipe for star formation is necessary. We deflne the stellar birthrate function as the 
number of stars formed in the time interval dt and in the mass range dm as: 

B(ra, t) = ■)Jj{t)Lp{m)dtdm (1.1) 

where: 

^j{t) = SFR (1.2) 

is the star formation rate (SFR), and: 

tp{m) = IMF (1.3) 

is the initial mass function (IMF) . The SFR is assumed to be only a function of time and 
the IMF only a function of mass. This is clearly an oversemplification but is necessary 
in absence of a clear knowledge of the star formation process. 

• Stellar evolution and nucleosynthesis- Nuclear burnings take place in the stellar 
interiors during the star lifetime and produce new chemical elements, in particular metals. 
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These metals, together with the pristine stellar material is restored into the interstellar 
medium (ISM) at the star death. This process clearly affects crucially the chemical 
evolution of the ISM. In order to take into account the elemental production by stars we 
define the "yields" , in particular the stellar yields (the amount of elements produced by 
a single star) and the yields per stellar generation (the amount of elements produced by 
an entire stellar generation). 

• Supplementary Parameters- Infall of cxtragalactic gas, radial flows and galactic 
winds are important ingredients in building galactic chemical evolution models. 



2. The stellar birthrate 

2.1. Theoretical recipes for the SFR 

Several paramctrizations, besides the simple one of a constant ilj{t), are used in the 
literature for the SFR: 

• Exponentially decreasing: 

SFR = ve-^l-"' (2.4) 

with T, =5—15 Gyr in order to produce realistic values which can be compared with 
the present time SFR in the Milky Way (Tosi, 1988). 

• The Schmidt (1965) law: 

SFR = vu\^, (2.5) 
is the most widely adopted formulation for the SFR. It was originally formulated by 
Schmidt (1959;1963) as a function of the volume gas density with k = 2.0. He measured 
the space density of stars in different regions of the Galaxy in relation to the number 
density of neutral hydrogen, measured by means of the 21 cm emission The formulation 
as a function of the surface gas density a gas is normally preferred for studying the Milky 
Way disk and galactic disks in general. In principle, Schmidt's formulation as a function 
of the volume gas density and that with a gas arc equivalent when fc = 1. Kennicutt, 
in a series of papers (1983;1989;1998a,b) tried to assess the dependence of massive star 
formation on the surface gas density in disk galaxies, by comparing Ha emission with the 
data on the distribution of HI and GO and found that the SFR can be well represented 
by a Schmidt law with k = 1.4 ± 0.15. The quantity v is the efficiency of star formation 
and is expressed in units of time~^. 

• A more complex formulation, which depends also upon the total surface mass density, 
was suggested by Dopita and Ryder (1994) and can be written as: 

SFR = va1:,,a^gi, (2.6) 

with cTtot being the total surface mass density. This kind of SFR is related to the feedback 
mechanism between the energy injected into the ISM by supernovae (SNe) and stellar 
winds and the local potential well. 

• In alternative to the Schmidt law, Kennicutt proposed also the following star for- 
mation law, which fits equally well the observational data: 

SFR = Qmi^gasCTgas OC R'^CTgas (2.7) 

with Q.gas being the angular rotation speed of the gas. 

2.2. The tracers of star formation 
The main tracers of star formation in galaxies arc: 

• Gounts of luminous supergiants in nearby galaxies under the assumption that their 
number is proportional to the SFR. 
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• The Ha and if/j flux from HII regions, which are ionized by young and hot stars, 
under the assumption that such flux is proportional to the SFR (Kennicutt 1998): 

SFR{MQyr-^) = 7.9 • IQ-^^Lff^ {ergs'^) (2.8) 

• From the integrated UBV colors and spectra of galaxies one can estimate the relative 
proportions of young and old stars and derive the ratio between the present time SFR 
and the average SFR in the past. 

• The frequency of type II SNe as well as the distribution of SN remnants and pulsars 
can be used as tracers of the SFR. These tracers have been used for deriving the SFR in 
the Galactic disk. 

• The radio emission from HII regions can also be a tracer of the SFR. 

• The ultraviolet continuum and the infrared continuum (star forming regions are 
surrounded by dust) are also connected to the SFR as in the following expression: 

SFR{MQyr-^) = 0.9 ■ 10"'^ '^^^^^ (2.9) 

EbolQ 

derived by means of a two-slope IMF by Donas et al. (1987). 

• Finally, the SFR can be derived from the distribution of molecular clouds (Rana 
1991). 

All of these formulations for the SFR need the assumption of an IMF and viceversa, 
the derivation of the IMF needs the assumption of a star formation history, as described 
in the next section. The local SFR, derived under the assumption of a particular IMF, 
suitable for the solar neighbourhood, gives the following range of values: 

SFR = 2 - lOMopc-^Gyr-i (2.10) 

(see Timmes et al. 1995) 

2.3. The IMF: Various Parametrizations 

The IMF is a probability distribution function and is normally approximated by a power 
law, namely: 



(p{m) = am-(i+^) (2.11) 

which is the number of stars with masses in the interval m, m+dm. The IMF can be 
one-slope (Salpotcr (1955) x=1.35) or multi-slope (Scalo 1986,1998; Kroupa et al. 1993). 
The IMF is usually normalized as: 

m(fi{m)dm = 1 (2-12) 

The IMF is derived locally from the present day mass function (PDMF) which in turn is 
obtained by counting the Main Sequence stars per interval of magnitude. Then the star 
counts are transformed into number of stars per pc^ and then a mass-luminosity relation 
is adopted to pass from the luminosity to the mass. 



2.4. Derivation of the IMF 

As already mentioned, the IMF is derived by the observed PDMF, which is the current 

mass distribution of Main-Sequence (MS) stars per unit area, n(m). 

For stars (0.1 < m/niQ < 1) with lifetimes Tm > to (with to being the Galactic age), 
n(m) can be written as: 
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n{m) = / (p{m)tl){t)dt (2.13) 
Jo 

where to is the Galactic lifetime. These stars are all still on the MS. If (p{m) is constant 
in time, as it is usually assumed, then: 

n{m) = ip{m) < ip > ta (2-14) 

where < if) > is, the average SFR in the past. 

For stars with Tm «tG {in> 2Mq), we see on the MS only those born after the time 
{t = tG-Tm)- The PDMF is therefore: 

n{m) = I ip{m)il){t)dt (2.15) 

JtG-Tm 

Again, if (p(m) is constant in time: 

n{m) = ip{m)ip{tG)Tm (2.16) 

under the assumption that tpita) = V'(*G — Tm), where tpito) is the SFR at the present 

time, to, and Tm is the lifetime of a star of mass m. 

The IMF, (p{m), in the mass interval 1 — 2 Mq depends on the ratio: 

6fc) = |M (2.1,) 

It has been shown by Scalo (1986) that a good fit between the two portions of the <^(m), 
namely below IMq and above 2Mq, requires: 

0.5 < fe(fG) < 1.5 (2.18) 

which means that the SFR in the local disk should have varied in time less than a factor 
of 2 during the whole disk lifetime. 

2.5. The Infall Rate: Various Parametrizations 

The presence of infall of extragalactic gas in the chemical evolution of galaxies is de- 
manded by the G-dwarf metallicity distribution in the solar vicinity and by the existence 
of high velocity clouds infalling towards the galactic disk. The origin of this infalling gas 
on the Galaxy is not yet entirely clear and measurements of the metallicity of such gas 
are necessary to decide whether it originates in the galactic disk (galactic fountain) or if 
it has an extragalactic origin. 

Several parametrizations for the infall rate have been used so far: 

• Constant in space and time 

• Variable in space and time such as: 

IR = >l(i?)e-*/^(^) (2.19) 

with t{R) is constant or varying along the disk and A{R) is derived by fitting the present 
time distribution of atotiR^tc)- 



3. Nucleosynthesis 

During the Big Bang the light elements (D, "^ffe, '^He and ^L«) were produced. On 
the other hand, all the elements heavier than "^Li, with the exception of Be and B, are 
produced inside stars. The light elements ^Li, Be and B are instead manufactured by 
spallation processes in the ISM due to the interaction between cosmic rays and interstellar 
atoms. 
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3.1. Nucleosynthesis in the Big Bang 

I give here a brief summary of the main steps in the Big Bang nucleosynthesis. 

When the temperature in the universe was T = lO^^K, only weak interactions causing 
conversions between protons and neutrons occurred, namely: 

p + e^n + u, (3.20) 

p + v^n + e'^. (3.21) 

The nucleosynthesis started when T = lO^K and lasted until T = lO^K. The first 
element to be formed was D (a nucleus composed by a neutron plus a proton) and 
subsequently ^He from the reaction: 

D+p^^ He + -i (3.22) 

followed by: 

^He + n^^ He + j (3.23) 

Then also very small fractions of "^Li (10~^ by mass) and ^Be (10~^^ by mass) were 
produced. 

One of the major achievements in cosmology is that it can account simultaneously 
for the primordial abundances of H, D, ^He, "^He and Li but only for a low density 
universe. The comparison between the observed primordial abundances and the Big 
Bang nucleosynthesis calculations can allow to impose constraints upon the baryon to 
photon ratio (77) in the imivcrsc. In particular, for a baryon to photon ratio ~ 3 • 10""^*^ 
the baryonic density parameter of the universe is (Peacock, 1999): 

0.010 < < 0.015 (3.24) 

3.2. Stellar Nucleosynthesis 

Before discussing stellar nucleosynthesis we need to define the crucial stellar mass ranges. 

• Brown Dwarfs (M < Ml, Ml = 0.08 - O.O9M0). They never ignite H and their 
lifetimes are larger than the age of the universe. 

• Low mass stars (0.5 < M/Mq < Mhcf) (MHeF=l-85-2.2M0 depending on stellar 
models) ignite He explosively and become C-0 white dwarfs (WD). If M < O.5M0 they 
become He WD. The lifetimes range from several 10® years up to several Hubble times. 

• Intermediate mass stars {MhbF < M/Mq < M^p) ignite He quiescently. M„p is 
the limiting mass for the formation of a C-0 degenerate core and is in the range S-QMq, 
depending on stellar evolution models. Their lifetimes range from several 10^ to 10® 
years. These stars die as C-0 WDs if they are not in binary systems. If in binary 
systems, stars in this mass range can give rise to Type la SNe (see later). 

• Massive stars {M > Mup) 

Stars in the mass range M^p < M/Mq < 10 — 12 become e-capture SNe (Type II SNe) 
and leave neutron stars as remnants. Stars in the range 10 — 12 < M/Mq < M^r 
{MwR ~ 20 — 40Mq) end their lives as core-collapse SNe (Type II) and leave a neutron 
star or a black hole as a remnant. Stars in the range Mwr < M/Mq < 100 probably 
become Type lb SNe. The lifetimes of these stars are in the range from several lO'' to 
^10^ years. 

• Very Massive Stars (M > lOOM©), if they exist, explode by means of "pair creation" 
and are called pair-creation SNe. In fact, at T ~ 2 • 10® K a large portion of the 
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gravitational energy goes into creation of pairs (e"*", e~), thus the star becomes unstable 
and explodes. These SNe leave no remnant and have lifetimes < 10^ years. 

• Supermassive objects (400 < M/Mq < 7.5 • 10^), if they exist, either explode due to 
explosive H-burning or collapse directly to black holes. The only available nucleosynthesis 
calculations for these stars are from Woosley ct al. (1984). 

All the elements with mass number A from 12 to 60 have been formed in stars during 
the quiescent burnings occurring during their lifetime. The main nuclear burnings are 
H, He, C, Nc, O and Si. Stars transform H into He and then He into heaviers until 
the Fe-peak elements, where the binding energy per nucleon reaches a maximum. At 
this point nuclear fusion reactions cannot occur anymore and the Fe nucleus starts con- 
tracting. When the central density reaches the atomic nuclear density, matter becomes 
uncompressible and a core-bouncx' occurs with the consequent formation of a shock wave 
and ejection of the star mantle. However, a problem exists for the explosion of stars with 
large Fe cores since most of the collapse gravitational energy is used to photodisintegrate 
Fe, thus weakening the shock wave and preventing the mantle ejection. To overcome 
this problem it has been suggested the existence of some mechanisms able to rejuvenate 
the shock wave, such as neutrino-heating from the collapsing neutron star, rotation and 
magnetic fields. 

Here is a summary of the main nucleosynthesis stages: the first element to be burned is 
H, which is transformed into He through the proton-proton chain or the CNO-cycle, then 
'^He is transformed into ^^C through the triple a reaction. Elements heavier than ^^C are 
then produced by synthesis of a-particles thus producing the so-called a-elements (O, Ne, 
Mg, Si, S, Ca and Ti). The last main burning in stars is the "^^Si -burning which produces 
^^Ni which then /3-decays into ^^Co and ^^Fe. Si-burning can be quiescent or explosive 
(depending on the temperature) but it always produces Fe. Explosive nucleosynthesis 
occurs in the inverse order (Si, O, Ne, C, Ho, H) relative to quiescent nucleosynthesis, 
depending on the fact that it starts in the center and propagates outwards following 
the passage of the shock wave. The main products of explosive nucleosynthesis are the 
Fe-peak elements. It is worth noting that low and intermediate mass stars never ignite 
C and thus end their lives as C-0 white dwarfs. Stars with masses below 10-12 Mq 
(depending on stellar models) explode during 0-burning and end up as e-capture SNe. 
Only massive stars can ignite all six nuclear fuels until they form an Fe-core. S- and 
r-process elements (elements with A> 60 up to Th and U) are formed by means of slow 
or rapid (relative to the f3- decay) neutron capture by Fe seed nuclei. In particular, s- 
processing occurs during quiescent He-burning both in massive and low and intermediate 
mass stars, whereas r-processing occurs during SN explosions. 

3.3. Supernova Progenitors 

Supcrnovac, planetary nebulae (PNc) and, to a minor extent, stellar winds arc the means 
to restore the nuclearly enriched material into the ISM, thus giving rise to the process of 
chemical evolution. There are two main Types of SNe (II, I) then divided in subclasses: 
SNe IIL, IIP and SNe la, lb, Ic. As already mentioned before, SNe II, which are believed 
to be the end state of stars more massive than IOMq exploding after a Fe core is formed 
(core-collapse SNe), produce mainly a-elements (O, Ne, Mg, Si, S, Ca) plus some Fe. 
The amount of Fe produced by type II SNe is one of the most uncertain quantities since 
it depends upon the so-called mass cut (how much Fe remains in the collapsing core and 
how much is ejected) and on explosive nucleosynthesis. 

Type la SNe arc believed to originate from the C-deflagration of a WD reaching the 
Chandrasekhar mass (1.44 Mq) after accretion of material from a young companion in 
a close binary system. C-deflagration occurs as a consequence of such accretion and 
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an explosion ensues destroying the whole star. No remnant is then left behind. They 
produce a large amount (~ 0.6 — O.TM©) of ^^Ni (Fe) plus traces of C to Si elements. 
C-deflagration is the explosive burning which best reproduces the observed abundance 
pattern in Type la SN remnants. The best model for this kind of explosion is model W7 
by Nomoto, Thielemann & Yokoi (1984). The amount of Fe produced by the other Type 
I SNe is smaller than that produced by the Type la ones at least a factor of two or more. 
Therefore, Type la SNe should be considered as the responsible for the Fe enrichment in 
the universe. 

3.4. Element production 

Here is a summary of element production: 

• Big Bang light elements H, D, ^He, ^i?e, '^Li. Deuterium is only destroyed inside 

stars to form "^He. ^He is also mainly destroyed. The only stars producing some '^He 
arc those with masses < 2.5Mq. Recent prescriptions for the yields of ^He arc from 
Forcstini & Charbonncl (1997) and Sackmann & Boothroyd (1999). Lithimii: ''Li is 
produced during the Big Bang but also in stars: massive AGB stars, SNe II, carbon- 
stars and novae. Some '^Li should also be produced in spallation processes by galactic 
cosmic rays (sec Romano et al. 2001). 

• Spallation Processes — > ^Li, Be and B. 

• Type II SNe a-elements (O, Ne, Mg, Si, S, Ca), some Fe, s-process elements {A < 
90) and r-process elements. Yields are from Woosley & Weaver (1995) and Thielemann 
et al. (1996). 

• Type la SNc Fc-peak elements. Yields arc from Nomoto ct al. (1984) and Nomoto 
et al. (1997). 

• Low and intermediate mass stars ^ *He. C. N, s-process {A > 90) elements. Yields 
are from Renzini & Voli (1981), Marigo et al. (1996), van den Hoek & Groenewegen 
(1997) and Gallino et al. (1998). 

3.5. Stellar yields 

In order to include the results from nucleosynthesis into chemical evolution models we 

need to define the stellar yields. The stellar yield of an element i is defined as the mass 
fraction of a star of mass m which has been newly created as species i and ejected: 

Pim = (^)i (3.25) 
m 

In order to compute pim we need to know some fundamental quantities from stellar 
evolution and nucleosynthesis: 

• Ma is the mass of the He-core (where H is turned into He) 

• Mco is the mass of the C-0 core (where the He is turned into heaviers) 

• Mrem is thc mass of the remnant (WD, neutron star, black hole) 
These masses are related to each other by thc following relations: 

Ma - Mhc = Mco 
where Mue is the newly formed and ejected '^He and: 

Mco -Mr^ Mc + Mo + Mheavters 

The values of these different quantities are given by stellar evolution and nucleosynthesis 
calculations. 
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4. Modelling chemical evolution 

4.1. Analytical models 

The simplest model of chemical evolution is the Simple Model for the chemical evolution 
of the solar neighbourhood. The basic assumptions of the Simple; Model are: 

- the system is one-zone and closed, namely there are no inflows or outflows, 

- the initial gas is primordial (no metals), 

- (/?(m) is constant in time, 

- the gas is well mixed at any time. 

In the following we will adopt the formalism of Tinsley (1980) and define: 

M = §^ (4.26) 

as the fractional mass of gas, with: 

Mtot = + Mgas (4.27) 

where M* is the mass in stars (dead and alive). Possible non-baryonic dark matter is 

not considered. 

The mass of stars can be expressed as: 

= (l-/x)Mtot. (4.28) 

The abundance by mass of an element i is defined by: 

X. = ^ (4.29) 

whore Mi is the mass in the form of the specific clement i. It is well known that the 
abundances must satisfy the condition, X)j^i(t) = 1, where the summation is over all 
the chemical elements. 
The initial conditions are: 

Mgas{^)=Mtot ;X,(0) = 0, (4.30) 

where i refers to metals. The equation for the evolution of the gas in the system can be 
written as: 

^E{t) (4.31) 

where E{t) is the rate at which dying stars restore both the enriched and unenriched 
material into the ISM. E{t) can be written as: 

/•oo 

E{t) = / (m - Mrem)^{t - T^)ip{m)dm (4.32) 

Jm{t) 

where m — Mrem is the total mass ejected from a star of mass m, and Tm is the lifetime 
of a star of mass m. When E{t) is substituted into the gas equation one obtains an 
integer-differential equation which can be solved analytically only by assuming Instan- 
taneous Recycling Approximation (I.R.A.). In this approximation one assumes that all 
stars less massive than IMq live forever whereas all stars more massive than IM© die 
instantaneously. In other words, I.R.A. allows us to neglect the stellar lifetimes and solve 
analytically equation (4.31). Under I.R.A. , we can define the returned fraction: 

/oo 
(to - Mrem)v{m)dm (4.33) 

which is called fraction because is divided by Jj°° mip{m)dm = 1, which depends on the 
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normalization of the IMF. We also define the yield per stellar generation as: 

1 f°° 

Vi = 5 / mpim(p{m)dm (4.34) 



where is the stellar yield previously defined. Then, by substituting R and t/j into 
equation (4.32) we obtain: 

E{t) = i){t)R (4.35) 

and: 

^ = -V(t)(l-i?) (4.36) 
The equation for the evolution of the chemical abundances can be written as: 

d{XiMgas) 



dt 

where: 



= -XiV(0 +Ei{t) (4.37) 



^ii't) = [{m- Mrem)Xi{t -Tm) +mpim]- tpit - Tm)ip{m)dm (4.38) 

Jm(t) 

contains both the unprocessed and the newly produced clement i. It is worth noting 
that this equation is valid for metals and not for elements which are wholly or partly 
destroyed in stars. Under the assumption of I.R.A. the above eq. becomes: 

Ei{t) = xP{t)RXi{t) + yi{l - R)ij{t) (4.39) 

When substituted in (4.37) the equation can be solved analytically with the previous 
initial conditions and the solution is: 

Xi = yM-) (4.40) 
A* 

the famous solution for the Simple Model. The yield which appears in the above solution 
is known as effective yield, simply defined as the yield j/i^j,^ that would be deduced if the 
system were assumed to be described by the Simple Model: 



y^.ff = TTTT^ (4-41) 



ln{l/y) 

The meaning of the effective yield can be understood with the following example: if 
Viefi > yi(true yield) then the actual system has attained a higher abundance for the 
element i at a given gas fraction /x. 

4.2. Failure of the Simple Model 

The Simple Model predicts too many stars with metallicity lower than [Fe/H]= -1.0 
dex relative to observations. This is known as "G- DWARF PROBLEM" . However, the 
G-dwarf is no more a problem since several solutions have been suggested. 

Possible solutions to the G-dwarf problem include: 

-Slow formation of the solar vicinity by gas infall 

-Variable IMF 

-Pre-enriched gas 

Generally, the slow infalling gas is preferred since it is the most realistic suggestion 
and produces results in very good agreement with the observations, as we will see in the 
following. 
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4.3. Analytical models with gas flows 
The equation for the evolution of abundances in presence of gas flows transforms into: 

d{XiMgas) ^ _x^(t)ij{t) + Ei{t) + XA,A{t) - Xi{t)W{t) (4.42) 

where A{t) is the accretion rate of matter with abundance of the element i X^i and W{t) 
is the rate of loss of material from the system. The case A{t) = W{t) = obviously cor- 
responds to the Simple Model. The case A{t) = 0, W{t) ^ corresponds to the outflow 
model. The easiest way of defining W{t) in order to solve the equation analytically is to 
assume: 

W{t) = A(l - R)ij{t) (4.43) 

where A > is the wind parameter. The analytical solution for the equation of metals, 
which can be integrated between e Xi{t) and between Mgas{0) = Mtot and Mgas{t), is: 

Xi = J^f'H^ + A)m-' - A] (4.44) 

For A = we recover the solution of the Simple Model. 

The case of A{t) ^ and W{t) = corresponds to the accretion model. The easiest 
way to choose the accretion rate is: 

A{t) = A(l - R)i:{t) (4.45) 

with A being a positive constant different from zero. The solution of the equation of 
metals for a primordial infalling material (X^i = 0) and A 7^ 1 is : 

= - (A - (A - l)/i-i)-^/(i-^)] (4.46) 

as shown by Matteucci and Chiosi (1983). If A = 1 ( extreme infall model) the solution 
is: 

X,=y,[l-e-(''"'-i)] (4.47) 

where the quantity /x"^ — 1 represents the ratio between the accreted mass and the initial 
mass. 



5. Equations with Type la and II SNe 

In general, if one wants to compute in detail the evolution of the abundances of elements 
produced and restored into the ISM on long timescales, the I.R.A. approximation is a bad 
approximation. Therefore, it is necessary to consider the stellar lifetimes in the chemical 
evolution equations and solve them with numerical methods. 

If Gi is the mass fraction of gas in the form of an element i, we can write: 

Gi{t) = -mMt) 

l-MBm 

+ / ll){t - Tm)Qmi{t - Tm)(t){m)dm 

J Ml 

pMbm 

-\-A I (t){m) 

JMBm 
/■0.5 

■[/ f{fJ')'^{t-Tm2)Qmiit-Tm2)dfj]dm 

/•Mbm 

+B / ■lp{t - Tm)Qmi{t - Tm)(t){m)dm 

•I Mbtt^ 
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+ I ll){t - Tm)Qmi{t - Tm)(l>{m)dm. 

J Mbm 

+XA,A{t) - XiW{t) (5.48) 

where B=l-A.The parameter A represents the unknown fraction of binary stars giving rise 
to type la SNe and is fixed by reproducing the observed present time SN la rate. Generally, 
values in the range A = 0.05 — 0.09 (according to the IMF) reproduce well the observed 
SN la rate in the Galaxy. The chemical abundances are defined as: Xj = where 

G = fi = = The total mass Mtot (or surface mass density atot) refers to the 

mass of stars (dead and alive) plus gas at the present time. 

Wi{t) = XiW{t) is the galactic wind rate for the element i, whereas Ai{t) = A{t)XAi, 
with A{t) being the accretion rate for the clement i. 

Mb m is the total minimum and Mbm the total maximum mass allowed for binary 
systems giving rise to Type la SNe (Matteucci & Greggio 1986). For the model of a C-0 
WD plus a red- giant companion for the progenitors of Type la SNe, Mbm < 16Mq . Mbwi 
is more uncertain and often has been taken to be ZMq to ensure that the primary star 
(the initially more massive in the binary system) would be massive enough to guarantee 
that after accretion from the companion the C-0 white dwarf eventually reaches the 
Chandrasekhar mass and ignites carbon. This formulation of the SN la rate was originally 
proposed by Greggio and Renzini (1983). Greggio (1996) presented revised criteria for 
the choice of Mbw.- In particular, the suggested condition for the explosion of the system 
is: 

MwD + eM2,e > Mch (5.49) 
(where M2,e is the envelope mass of the evolving secondary, Mwd is the mass of the 
white dwarf and Mch is the Chandrasekhar mass). 

For models involving Sub-Chandrasekhar white dwarf masses, which have been sug- 
gested to explain subluminous Type la SNe Greggio obtains: Mwd > O.6M0 and 
eM2,e > 0.15 

The masses Ml = 0.8 and Mu = lOOM© define the lowest and the highest mass, 
respectively, contributing to the chemical enrichment. The function Tm(rn) describes the 
stellar lifetimes. The quantity Qmi{t — t„i) contains all the information about stellar 
nucleosynthesis for elements either produced or destroyed inside stars or both, and is 
defined as in Talbot and Arnett (1973). 

5.1. Type la SN rates 

The single degenerate (SD) scenario is based on the original suggestion of Whelan and 
Iben (1973), namely C-deflagration in a WD reaching the Chandrasekhar mass after 
accreting material from a star which becomes red giaut and fills its Roche lobe. An 
alternative to the SD scenario is represented by the double degenerate (DD) scenario, 
where the merging of two C-0 WDs, due to gravitational wave radiation, creates an object 
exceeding the Chandrasekhar mass and exploding by C-defiagration (Iben and Tutukov 
1984). Negative results from observational searches for very close binary systems made 
of massive enough WDs (Bragaglia et al. 1990) has made the DD scenario less attractive 
and people to concentrate more on the SD scenario. 

A recent model has been suggested by Hachisu et al.(1996; 1999) and is based on the 
classical scenario of Whelan and Iben (1973) but with a metallicity efi'ect. It predicts that 
no Type la systems can form for [Fe/H]< —1.0. This model seems to have some difficulty 
in explaining the low [a/Fe] ratios observed in Damped Lyman- a systems (DLAs) which 
show that even at low metallicities is present the effect of Type la SNe. Recently, 
Matteucci & Recchi (2001) showed that there are some problems with this scenario also 
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in explaining the observed [0/Fe] vs. [Fe/H] relation in the solar neighbourhood and 
concluded that the best model is still the SD one in the formulation proposed by Greggio 
& Renzini (1983). They also showed that the typical timescale for Type la SN enrichment, 
defined as the time tsNia when the SN rate reaches the maximum, varies strongly from 
galaxy to galaxy and that it is not correct to adopt a universal tsNia=^ Gyr, as often 
quoted in the literature. Their results indicate that for an elliptical galaxy with high 
SFR, tsNia = 0.3 - 0.5 Gyr, for a spiral Galaxy like the Milky Way, tsNia = 4-5 Gyr 
and for an irregular galaxy with a continuous but very low SFR, tsNia > 5 Gyr. This 
fact has very important consequences on the chemical evolution of galaxies of different 
morphological types, as we will see in the following. 

6. The formation and evolution of the Milky Way 

The first application of models of chemical evolution is the Milky Way (MW) Galaxy 
for whic;li we have most of the available data. Before describing the chemical evolution 
of the Galaxy 1 recall some of the ideas proposed for the formation of the MW. Eggen, 
Lynden-Bell & Sandage (1962) suggested a rapid collapse for the formation of the Galaxy 
lasting ^ 3 • 10* years implying that no spread in the age of Globular Clusters should 
be observed. They based their suggestion on the finding that halo stars have high radial 
velocities whitnessing the initial fast collapse. Later on, Searle & Zinn (1978) proposed a 
central collapse but also that the outer halo formed by merging of large fragments taking 
place over a considerable timescale > 1 Gyr. More recently, Berman & Suchov (1991) 
proposed the hot Galaxy picture, an initial strong burst of star formation which inhibited 
further star formation for a few gigayears, while a strong galactic wind was created. 
Subsequently, the remainder of the proto-Galaxy, contracted and cooled to form the 
major stellar components observed today. At the present time, the most popular idea 
on the formation of the Galaxy is that the inner halo formed rather quickly on a time 
scale of 0.5-1 Gyr, whereas the outer halo formed more slowly by mergers of fragments 
or accretion from satellites of the MW (see Matteucci 2001). 

6.1. Models for the Milky Way 

• SERIAL FORMATION APPROACH: halo, thick and thin disk formed in sequence, 
as a continuous process (e.g. Matteucci & Frangois 1989). 

• PARALLEL FORMATION APPROACH: the various Galactic components start 
forming at the same time and from the same gas but evolve at different rates (e.g. Pardi, 
Ferrini & Matteucci 1995). At variance with the previous scenario it predicts overlapping 
of stars belonging to the different components. 

• TWO-INFALL APPROACH: the evolution of the halo and disk are totally indepen- 
dent and they form out of two separate infall episodes (overlapping in metallicity is also 
predicted) (e.g. Chiappini, Matteucci & Gratton, 1997; Chang et al. 1999, Alibes et al. 
2001). 

• STOCHASTIC APPROACH: in the early halo phases, mixing was not efficient, thus 
pollution from single SNe (Tsujimoto et al. 1999; Argast et al. 2000; Oey 2000) can be 
seen in very metal poor stars. This approach predicts a large spread in the abundance 
ratios at very low [Fe/H], even larger than observed. 

6.2. The two-infall model 
Here I describe in more detail the two-infall approach (Chiappini et al. 1997) since it 
gives the best agreement with observations for the halo and disk. 

The basic equations for the evolution of the MW in this case contain only an accretion 
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These equations are the same as (5.48) plus a term which contains the contribution 
from novae, which are perhaps important producers of "^Li, ^^N and ^^G (D'Antona 
& Matteucci, 1991;Romano et al. 2001). The nova contribution contains the parameter 
a=0.0155 which represents the fraction of WDs which are in binary systems giving rise to 
nova events. The value of a is chosen to reproduce the present time nova rate Rnovae ~ 
26yr~^, after assuming that each nova has ~ 10^ outbursts all over its lifetime. The 
quantity At ^ 1 Gyr is the assumed time-delay for the starting of nova activity since the 
formation of the binary system. 

The infall rate term is defined as: 



^ dGi{r,t) ^ A{r){X^)ie-'/-- ^ B(r)(X^),e-(*-*— )/-°('-) 
dt (Ttot{r,tG) otot{r,tG) 

where tr is the timescale for the inner halo formation (0.5-1 Gyr) and Toir) is the 
thin disk timescale varying with galacto centric distance (inside-out formation): 



Tr,(r) = 0.875r- 0.75 (6.51) 

The SFR is given by: 

m=vatot{r,tp<Jaas{tf^ (6.52) 
with k\ =0.5 and k2 = 1.5. The behaviour of this SFR for the halo- thick disk and the 

thin-disk phase, respectively, is show in Figure 1. A threshold density {ath = 7MqPC~'^) 
in the SFR is assumed. 

6.3. Applications to the Local Disk 

• The G-dwarf metallicity distribution is shown in Figure 2 where the data are com- 
pared with the predictions of the two infall model assuming a time scale for the formation 
of the local disk of 8 Gyr. (Chiappini et al. 1997; Boissier and Prantzos 1999; Chang et 
al. 1999; Chiappini et al. 2001). 

• The relative abundance ratios as functions of the relative metallicity (relative to the 
Sun) fX/Fe] vs. [Fe/H], they arc interpreted as due to the time-delay between Type la 
and 11 SNe. In fact, for the a-elements the slowly declining [a/Fe] ratio at low [ Fe/H] 
(halo phase) is normally interpreted as due to the pollution from massive stars, whereas 
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Figure 1 . Predicted star formation rate in the halo-thick disk phase and in the thin disk phase 
during the evolution of the Milky Way. Notice that because of the existence of a threshold in 
the gas density the SFR halts in between the two major episodes of infall and oscillates in the 
last phases of the evolution of the disk. The models are from Chiappini et al. (2001). 



the abrupt change in slope occurring at ~ \Fe/H] = — 1.0 dex is due to the Type la SNe 
restoring the bulk of iron. From the [a/Fe] vs. [Fe/H] diagram one can infer the timescale 
for the formation of the halo (r/j ^ 1.5-2.0 Gyr, Matteucci and Frangois, 1989; Chiappini 
et al. 1997), just by means of the age- [Fe/H] relationship, which indicates the time at 
which the metallicity of the turning point is reached. In Figure 3 we do not show the 
usual plot but [Fe/0] vs. [0/H] since in this plot are evident some features which do not 
appear in the classical diagram. In particular, the data show evidence for a gap around 
[0/H]== -0.3 dex, corresponding to [Fe/H] ~ -1.0 dex (Gratton et al. 2000). This gap is 
well reproduced by the model which predicts a hiatus (no more than 1 Gyr) in the SFR 
between the end of the thick-disk phase and the beginning of the formation of the thin 
disk. This hiatus produces, in fact, a situation where O is no longer produced whereas Fe 
is produced, and this is revealed by the increase of [Fe/a] at constant [a/H]. This effect 
has been observed also for [Fe/Mg] vs. [Mg/H] by Furhmann (1998). The model shown 
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Figure 2. Observed and predicted G-dwarf motallicity distribution. The data are from 
Rocha-Pinto and MacicI (1996) and Hou et al. (1998). The model predictions (continuous 
hue) are from Chiappini et aL (2001). The model assumes that the timescale for the disk 
formation in the solax neighbourhood is 8 Gyr. 



in Figure 3 is a two-infall model with a threshold density for the star formation and is 
just the existence of such a threshold which produces the hiatus in the SFR, evident in 
Figure 1. 

6.4. Applications to the whole disk 

• Abundance Gradients are known to exist along the Galactic disk from data from 
various sources (HII regions, PNc, B stars). These data suggest that the gradient for 
oxygen is ~ —0.07 dex/kpc in the galoctocentric distance range 4-14 kpc. It is not yet 
clear if the slope is unique or if there is a change in slope as a function of the galactocentric 
distance. Similar gradients are found for N and Fe (see Matteucci 2001 and references 
therein) . 

• Gas Distribution. HI is roughly constant over a range of 4-10 kpc along the Galactic 
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Figure 3. The [Fo/0] vs. [O/H] relation in the solar neighbourhood. The data are from 
Gratton et al.(2000) and the model (contirmous line) from Chiappini et al.(2001). In this figure 
is evident the existence of a gap in the data at around [O/H]=-0.3 dex, corresponding to [Fe/H] 
-1.0 dex, which is also predicted by the model. 



disk -while H2 follows the light distribution. No models can explain the two distributions. 

The total gas increases towards the center with a peak at 4-6 kpc. 

• The SFR Distribution is obtained from various tracers (Lyman-a continuum, pul- 
sars, SN remnants, molecular clouds) and shows that the SFR increases with decreasing 
galactoccntric distance reaching a peak at 4-6 kpc in correspondence of the gas peak. 

In order to fit gradients, SFR and gas one has to assume that the disk formed inside- 
out, in agreement with a previous suggestion by Larson (1976) and that the SFR should 
be a strongly varying function of the galactoccntric distance. In Figure 4 is shown a 
comparison between an inside-out model and the abundance gradient, the gas, star and 
SFR distributions. 
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Figure 4. Comparison between data and model predictions from Chiappini et al. (2001) 
(continuous lines). In the first panel from the top we show the oxygen abundance gradient. In 
the second panel the variation of the SFR/ SFRq, in the third panel the gas distribution and 
in the fourth panel the distribution of stars along the Galactic disk. 

6.5. The Role of Radial Flows in the evolution of the Galactic Disk 

The gas infalling onto the disk can induce radial inflows by transferring angular nionien- 
tum to the gas in the disk. Angular momentum transfer can be due to the gas viscosity 
in the disk and induce inflows in the inner parts of the disk and outflows in the outer 
parts. All viscous models suggest that mctallicity gradients can be steepened by radial 
(in)flows, especially if an outer star formation cut-off is assumed (Clarke 1989; Yoshii & 
Sommer-Larsen 1989). All models agree that the velocity of radial (in)flows should be 
low (f < 2 km/sec). Obscrvationally is not clear if these radial flows exist. Edmunds & 
Greenhow (1995), by means of analytical models of galactic chemical evolution concluded 
that there is no simple one-way effect of radial flows on abundance gradients. Portinari 
& Chiosi (2000) suggested that radial inflows may represent a possible explanation of 
the peak of the gas at 4-6 kpc. In our opinion radial flows, if they exist, are never the 
main cause for the formation of abundance gradients. 
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6.6. The Role of the IMF in the evolution of the Galactic Disk 

Abundance gradients, in principle can be obtained by a variation of the IMF along the 
disk. One can assume either that more massive stars form in external regions or the 
contrary, that more low mass stars form in external regions. It has been shown convinc- 
ingly that neither of the two options work (Carigi 1996; Chiappini et al. 2000), and that 
there is no evidence in favor of such variations along the Galactic disk. Therefore, we can 
exclude the variation of the IMF as the main cause of abundance gradients and conclude 
that the best agreement with the observed properties of the Galactic disk is obtained by 
assuming a constant IMF. 

6.7. Scenarios for Bulge Formation 

Various scenarios have been proposed insofar for the formation of the Galactic bulge but 
only one seems to reproduce the observed abundance pattern. Here I recall the main 

scenarios: 

• Accretion of extant stellar systems which eventually settle in the center of the 
Galaxy. 

• Accumulation of gas at the center of the Galaxy and subsequent evolution with 
either fast or slow star formation. 

• Accumulation either rapid or slow of metal enriched gas from the halo or thick disk 
in the Galaxy center. 

• Formation occurs out of inflow of metal enriched gas from the thin-disk. 

The metallicity distribution of stars in the bulge as well as the [a/Fe] vs. [Fe/H] relations 
help in selecting the most probable scenario and suggest that a fast accumulation of gas 
in the Galactic center accompanied by fast star formation is the best scenario. In this 
framework, the bulge must have formed contemporarily to the inner halo on a similar 
timescale. In Figures 5 and 6 we show some model predictions compared with the 
metallicity distribution of bulge stars and the predicted [a/Fe] vs. [Fe/H] ratios for the 
bulge together with the predictions for the solar neighbourhood. The model for the 
bulge presented in the two figures (Matteucci, Romano & Molaro, 1999) assumes a much 
stronger star formation rate than in the solar neighbourhood (by a factor of 10) with 
the same nucleosynthesis prescriptions and a timescale for bulge formation of 0,5 Gyr 
as opposed to 8 Gyr in the solar vicinity. In Figure 5 it is evident that the best model 
to reproduce the stellar metallicity distribution requires a Salpeter (1955) IMF which is 
sligthly flatter than that used for the solar neighbourhood (Scalo, 1986). The predicted 
[a/Fe] ratios in Figure 6 indicate that the bulge stars should show overabundances of 
a-elements for most of the [Fc/H] range, as it seems also suggested by observations (Mc 
William & Rich, 1994; Barbuy, 1999). This is a consequence of the time-delay between 
Type la and II SNe, coupled with a very fast evolution in the bulge as compared to the 
solar vicinity. In fact, in this case high values of [Fe/H] are reached in the gas before 
a substantial number of SNe la has the time to restore the bulk of iron. The contrary 
occurs in a system with lower star formation than in the solar neighbourhood, such as 
the external regions of the disk and irregular galaxies. In these systems we expect that 
the overabundance of a-elements relative to Fe is maintained only for a short interval of 
[Fe/H] (see Pagel, 1997: Matteucci 2001). 

7. Disks of Other Spirals 

Abundance gradients 

Abundance gradients in dex/kpc are known to exist also in the disk of other spirals 
showing that they are steeper in smaller disks, but this correlation disappears if the 




Figure 5. Comparison between data (dotted lines) and model predictions for the metallicity 
distribution of bulge stars. The data are from Mc William & Rich (1994) and the models (con- 
tinuous lines) from Matteucci et al. (1999). Each panel corresponds to a model with a different 
IMF slope, as indicated. 

gradients are expressed in units of dex/scalelength, thus indicating the existence of a 
universal slope per unit scalclcngth (e.g. Garnctt, 1998). Another remarkable charac- 
teristic about abundance gradients in other spirals is that they appear to be flatter in 
galaxies with central bars, suggesting that the dynamical effect of the bar can influence 
the evolution of the disk. 
The SFR 

The SFR is measured mainly from Ha emission (Kennicutt, 1998) and implies a cor- 
relation with the total surface gas density {lil+H2) as discussed in section 2.1. 

Gas distributions 

Gas distributions, especially the HI distribution is known for a fair sample of spirals. 
Thc;re are indications of diff'erences between fleld and cluster spirals (e.g. Skillman et al. 
1996). 

Integrated colors 
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Figure 6. Predicted [a/Fc] vs. [Fe/H] relations for the bulge (continuous lines) and for the 
solar neighbourhood (dotted lines). The models are from Matteucci et al. (1999). 

Studies of integrated colors of spiral disks reflect the observed abundance gradients 
and are well reproduced by an inside-out formation of disks similar to what assumed for 
the Milky Way (Josey & Arimoto 1992; Jimenez et al. 1998; Prantzos & Boissier 2000). 

8. Conclusions on the Milky Way and other spirals 

The comparison between obscirvations and models for the Milky Way suggests that: 

• The disk of the Galaxy formed mostly by infall of primordial or very metal poor 
gas accumulating faster in the inner than in the outer regions (inside-out scenario, i.e. 
td{R) decreases with decreasing R). 

• In the framework of the inside-out scenario, the SFR should be a strongly varying 
function of the galacto centric distance. In particular, it should either depend on the total 
surface mass density (feed-back mechanism) or on the angular circular velocity of gas. 
Both formulations for the SFR are supported by observations. Under the assumption of 
a strongly varying SFR, the observed abundance gradients can be nicely reproduced. 
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• Radial flows probably are not the main cause of gradients but can help in reproducing 
the gas profile. 

• A constant IMF (in time and space) should be preferred, since variable IMFs have 
been tested and they cannot reproduce all of the observational constraints of the solar 
neighbourhood and the whole disk at the same time. 

• The properties of the disks of other spirals also indicate an inside-out disk formation. 

• In the framework of semi-analytical models of galaxy formation (Mo, Mao & White, 
1998), the evolution of galaxy disks can be described by means of scaling laws calibrated 
on the Galaxy with Vc and A as parameters (Jimenez et al. 1998; Prantzos & Boissier, 
2000), where Vc is a measure of the mass of the dark halo and A is a measure of the 
specific angular momentum of the halo. 

• Dynamical processes such as the formation of a central bar can influence the evolu- 
tion of disks and deserve more attention in the future. 

• Abundance ratios (e.g. [a/Fe]) in stars al largo galactocentric distances can give 
us a clue to interpret the formation of the disk and the halo (inside-out or outside- in). 
In fact, in a typical inside-out scenario we predict that the [a/Fe] ratios would decrease 
with galactocentric distance, due to the weaker and sometimes intermittent SF regime, 
whereas in an outside-in scenario wc would expect the contrary due to the faster evolution 
of the more external region which quickly consume the gas before Type la SNe have time 
to restore the bulk of Fe. 

• The predictions of chemical evolution models can be tested in a cosmological context 
to study the galaxy surface brightness and size evolution as a function of redshift. Roche 
et al. (1998) have already done that and suggested that a size and luminosity evolution, 
as suggested by the inside-out scenario fits better the observations. 

• Inside-out formation of the Galaxy is suggested also by the fact that the globular 
clusters of the inner halo are coeval (At ~ 0.5 Gyr, Rosenberg et al. 1999), whereas the 
age difference seems to increase in the outer halo. 

• The metallicity distribution of stars in the Bulge as well as the observed [a/Fe] ratios 
suggest a very fast formation of the Bulge (~ 0.5-0.8 Gyr), due a the very fast SFR, as 
predicted by succesfull chemical evolution models and by Elmegreen (1999). 



9. Elliptical Galaxies 

Elliptical galaxies are mainly found in galaxy clusters, they show a large range of 
luminosities and masses and are made by stars as old as those in globular clusters. No 
cold gas is observed in these galaxies but hot X-ray gas halos are present. Because of 
their large masses and metal content they are probably the most important factories of 
metals in the universe. 

9.1. Observational properties 

Here I recall the main observational features of elliptical galaxies: 

• The existence of Color-Magnitude and color - velocity dispersion (cTq) relations (col- 
ors become redder with increasing luminosity and mass; e.g. Bower et al. 1992) is inter- 
preted as a metallicity effect, namely as the fact that the metallicity decreases outwards 
similarly to what happens in spiral disks. 

• The existence of a Mg2-(Jo (where Mg2 is a metallicity index) relation reinforces the 
previous point (Bender et al. 1993; Bernardi et al. 1998; CoUess et al. 1999; Kuntschner 
et al. 2001). 

• Inside ellipticals, Mg2 correlates also with the escape velocity first shown 
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by Franx & Illingworth (1990), indicating that the magnesium index is larger where the 
escape velocity is larger. 

• M/ Lb increases by a factor of ~ 3 from faint to bright ellipticals implying a tilt of 
the fundamental plane of ellipticals (Bender et al. 1992). The fundamental plane is the 
particular plane occupied by these galaxies in the space defined by the stellar velocity 
dispersion, the effective radius and the surface brigthness. 

• Abundance gradients inside ellipticals have been measured by means of metallicity 
indices such as Mg2 and < Fe > (CaroUo et al. 1993; Davies et al. 1993; Kobayashi & 
Arimoto 1999). These gradients correspond roughly to a gradient in [Fe/H] of the stellar 
component of A[Fe/H]/Alogr ~ -0.3. 

The average metallicity of the stellar component in ellipticals is < [Fe/H] >*~ — 0.3dex 
(from -0.8 to +0.3). 

• By comparing synthetic metallicity indices with the observed ones some authors 
(Worthey et al. 1992; Weiss et al. 1995; Kuntschner et al. 2001) have suggested that 

the average stellar < [Mg/Fe] >* is larger than zero (from 0.05 to + 0.3 dex) in nuclei 
of giant ellipticals. Moreover, there is indication that < [Mg/Fe] >* increases with 
increasing CTq [Mgai) and luminosity (Worthey et al. 1992; Jorgensen 1999; Kuntschner 
et al. 2001). 

In particular the relation found by Kuntschner et al. (2001) is: [Mg/Fe]=0.30(±0.06)Zofifcro— 
0.52(±0.15). 

9.2. Form,ation of Ellipticals 

Several mechanisms have been suggested for the formation and evolution of elliptical 
galaxies, in particular: 

• Early monolithic collapse of a gas cloud or early merging of lumps of gas where 
dissipation plays a fundamental role (Larson 1974; Arimoto & Yoshii 1987; Matteucci & 
Tornambe 1987 ). In this scenario the star formation stops soon after a galactic wind 
develops and the galaxy evolves passively since then. 

• Bursts of star formation in merging subsystems made of gas (Tinsley & Larson 1979). 
In this picture star formation stops after the last burst and gas is lost via stripping or 
wind. 

• Early merging of lumps containing gas and stars in which some dissipation is present 
(Bender et al. 1993). 

• Merging of early formed stellar systems in a wide rcdshift range and preferentially 
at late epochs (Kauffmann et al. 1993). A burst of star formation can occur during the 
major merging where ~ 30% of the stars can be formed (Kauffmann 1996). 

The main difference between the monolithic collapse scenario and the hierarchical 
merging relies in the time of galaxy formation, occurring quite early in the former scenario 
and continuously in the latter scenario. As we will see, there are arguments either in 
favour of the monolithic or the hierarchical scenario. 

9.3. Formation of Ellipticals at low z 

Here I recall some of the main arguments in favor of the formation of ellipticals at low 

redshifts: 

• Relative large values of the index measured in a sample of nearby ellipticals 
which could indicate prolonged star formation activity up to 2 Gyr ago (Gonzalez 1993; 
Trager et al. 1998). 

• The tight relations in the fundamental plane are due to a conspiracy of age and 
metallicity in the sense that it should exist an age-metallicity anticorrelation implying 
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that the more metal rich galaxies are also younger (Ferreras et al. 1999; Trager et al. 
2000). 

• The apparent paucity of high luminosity ellipticals at z ~ 1 compared to now claimed 
by a series of authors (Kauffmann et al. 1993; Zepf, 1997; Menanteau et al. 1999). 

9.4. Formation of Ellipticals at high z 

Here I recall the arguments in favor of a formation of ellipticals at high redshift: 

• The tightness of the color-central velocity dispersion relation found for Virgo and 
Coma galaxies (Bower et al. 1992). If the formation of ellipticals were a continuous 
process we should expect a much larger spread in the galaxy colors for a given central 
velocity dispersion. In particular, the argument goes like that: from the observed color 
scatter one can derive t}j — tp ^ 2Gyr (where t}j is the Hubble time and tp the time of 
galaxy formation). litH = 15 Gyr then the youngest ellipticals must have formed ~ 13 
Gyr ago at z > 2 (Renzini, 1994). 

• The thinness of the fundamental plane for ellipticals in the same two clusters, in 
particular the M/L vs. M relation (Renzini & Ciotti 1993). 

• The tightness of the color-magnitude relation for ellipticals in clusters up to 2; ~ 1 
(Kodama et al. 1998; Stanford et al. 1998) 

• The modest passive evolution measured for cluster ellipticals at intermediate redshift 
(van Dokkum & Franx 1996; Bender et al. 1996). 

• Lyman-break galaxies at z > 3 where the SFR~ 50 — lQQMQyr~^ could be the 
young ellipticals (Steidel et al. 1996; 1998). 

• The strongly evolving population of Luminous Infrared Galaxies suggesting that 
they are progenitors of massive spheroidals (Blain et al. 1999; Elbaz et al. 1999). 

9.5. Models for ellipticals based on galactic winds 
Monolithic models assume that ellipticals suffer a strong star formation and quickly 
produce galactic winds when the energy from SNe injected into the ISM equates the 
potential energy of the gas. Then, star formation is assumed to halt after the development 
of a galactic wind. In this framework, the evolution of ellipticals crucially depends on 
the time at which a galactic wind occurs, taw- For this reason, it is extremely important 
to understand the SN feedback and star formation process. The condition for the onset 
of a wind can be written as: 

{Eth)lSM > Esgas (9.53) 
The thermal energy of gas due to SNe and stellar wind heating is: 

{Eth)iSM = EthsM + Eth^ (9.54) 

with 

EthsM = I esNRsN{t)dt (9.55) 
Jo 

and 

rt /-lOO 

Eth^ = / ip{m)ilj{t)eu,dmdt (9.56) 
Jo J12 

for the contribution from SNe and stellar winds, respectively. The quantity Rsn repre- 
sents the SN rate (II and la). The quantities esN = Vsn^o with Co = 10^^ erg where €„ 
is the typical SN blast wave energy , and = rj^E^ with = 10^^ erg (typical energy 
injected by a 2QMq star taken as representative), arc the efficiencies for the energy trans- 
fer from SNe II and la into the ISM. These efficiencies, r]w and r]sN, can be assumed as 
free parameters or be calculated from the results of the evolution of a SN remnant in the 
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ISM and the evolution of stellar winds. The SN feedback is, in fact, a crucial parameter. 
The formulation of Cox (1972) for the efhciency of energy injection from SNe derives 
from following the evolution of the shock wave produced by the explosion into an ISM 
with constant density: 

esN = 0.72eo erg (9.57) 

for tsN <tc = bJlO'^et^^^rio^^^^ years, where tc is the cooling time, and tsN is the time 
elapsed from the SN explosion, and: 

esN = 2.2eo{tsN/te)-°-''^ erg (9.58) 

for tsN > tc- 

With these prescriptions only few % of arc deposited into the ISM (sec Bradamante 
et al. 1998). However, multiple SN explosion should change the situation. Unfortunately, 
very few calculations of this type are available. An important point to consider is also 
that SNc la, which explode after Type II SNc, should provide more energy into the ISM 
than Type II SNe since they explode in an already formed cavity (see Recchi, Matteucci 
& D'Ercole, 2001). 

The total mass of the galaxy is expressed as Mtot{t) = M^,{t) + Mgas{t) + Mj_ark{t) 
with M^it) = M^{t) + Mgas{t) and the binding energy of gas is: 

Esgasit) = WL{t) + WLoit) (9.59) 

and: 

WUt) = _o.5G«l^ (9.60) 
rL 

represents the potential well due to the luminous matter, whereas: 

WLoit) = -GWLD^^^^^^^^^^ (9.61) 

is the potential well due to the interaction between dark and luminous matter, where 
wld ^ ^S{1 + 1.375") with S = rL/rr,(Bertin et al. 1992). 
The SFR is usually assumed to be: 

SFR = vMgas (9.62) 

where v oc (7 =-0.11, Arimoto & Yoshii 1987), owing to the fact that the star 

formation efficiency is just the inverse of the timescale for star formation: 

V = Tgl (9.63) 

and that: 

TSF « Tcoll OC Tff (9.64) 

with Tcoll and Tff being the collapse and the free-fall timescale, respectively. Since 
dynamical timescales are longer for more massive galaxies, the efficiency of star formation 
should decrease with galactic mass. The efficiency v coupled with the increase of the 
potential well as increases leads to the fact that more massive galaxies form stars for 
a longer period before suffering a galactic wind. This fact has been invoked for explaining 
the observed mass-mctallicity relation (Larson, 1974). However, this is at variance with 
the observed [Mg/Fe] vs. CTo relation, since in this scenario the more massive ellipticals 
should show the lowest [Mg/Fe] . 
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9.6. Failure of Larson's Model 



There are different ways of obtaining that [Mg/Fe] in the steUar component increases 
when the galactic mass Ml increases. These ways are: 

• DifFcrcnt timcscalcs for star formation (Worthcy ct al. 1992) in the sense that star 
formation should be more efficient in more massive galaxies {i/ oc M^). In this case, the 
situation could be such that a galactic wind occurs earlier in more massive systems, the 
inverse wind scenario (Matteucci 1994). 

• A variable IMF from galaxy to galaxy favoring more massive stars (Mg producers) 
in more massive galaxies (Worthey et al. 1992; Matteucci, 1994). 

• Different amounts and/or concentrations of dark matter as functions of Ml In par- 
ticular, less dark matter should be present in the most massive systems (Matteucci, 
Ponzone & Gibson, 1998). As a consequence of this, again galactic winds occur earlier 
in more massive objects. 

• A selective loss of metals: more massive systems loose more Fe relative to a-elements 
than less massive ones (Worthey et al. 1992). 

As mentioned before, Matteucci (1994) proposed a model that she called inverse wind 
model, where the efficiency of star formation is an increasing function of galactic mass, 
thus implying a shorter period of star formation in massive ellipticals. In fact, the 
efficiency of star formation is chosen in such a way that in massive ellipticals the galactic 
wind occurs before than in less massive ones. This produces the increase of the [Mg/Fe] 
ratio as a function of galactic mass. This approach bears a resemblance with the merging 
scenario of Tinsley & Larson (1979) where the efficiency of star formation was assumed 
to increase with the total mass of the system. In the inverse wind model a very massive 
elliptical of Ml = lO^^M© starts developing a wind before 1 Gyr from the beginning of 
star formation, whereas a small ellipticals form stars for a longer period. As a result, 
the average < [Mg/Fe] >* is larger in massive than in small ellipticals. The same effect 
can be obtained by varying the IMF in such a way that more massive ellipticals tend 
to form more; massive! stars. It is worth noting that this second hypothesis could also 
explain the tilt of the fundamental plane in M/Lb- On the other hand, the hypothesis 
of the variable dark matter does not produce relevant effects, as shown by Matteucci et 
al. (1998). The only assmnption which has not been tested quantitatively is the selective 
loss of metals. The main problem with all of these alternatives is, in any case, the lack of 
a good physical justification depending on the poorly known processes of star formation 
and SN feeback, and in the future some effort should be devoted to these fields. 

It is worth noting that the hierarchical clustering scenario for galaxy formation cannot 
produce a solution to the observed [Mg/Fe] trend in ellipticals. In fact, it rather predicts 
the contrary, as shown by Thomas (1999) and Thomas et al. (2002). In fact, in the hier- 
achical clustering the period of star formation in the most massive ellipticals is predicted 
to be the longest thus favoring low [Mg/Fe] ratios in massive objects, as shown in Figure 
7. 



The metallicity in elliptical galaxies always refers to the metal content of the stars, in 
particular of the stellar population dominating in the visual light. For this reason we 
should define the average stellar metallicity. In particular, the average metallicity of a 
composite stellar population averaged on the mass is: 



9.7. Averaged Stellar Metallicities 




Si 



(9.65) 



26 Francesca Matteucci: Chemical Evolution of Galaxies and Intracluster Medium 

0.4 



0.3 



ID 

^ 0.2 



0.1 







1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 



T 



• • • 




Hierarchical 
Models 



_l I I I I I I I I I I I I I I I I I L_ 







100 200 300 

G (kra/s) 



400 



Figure 7. Predicted [a/Fe] vs. velocity dispersion for elliptical galaxies. Data are compared 
with predictions obtained by adopting the star formation history assumed in hierarchical clus- 
tering models for galaxy formation. The models and the figure are from Thomas et al. (2002). 

where 1 refers to the specific time ti and Si is the total mass of stars ever born. This 
is the real average metallicity but in order to compare models and observations it is 
more appropriate to define the metallicity averaged on the visual light. In particular, the 
average metallicity of a composite stellar population averaged on the light is: 



(9.66) 



where Uij is the number of stars in the abundance interval Xi and luminosity interval 
Lvj. 

It is worth noting that < Xi >m is larger than < Xi >l for galaxies with Ml < 10^ Mq, 
since metal poor giants dominate the visual light whereas they are similar for larger 
masses (Yoshii & Arimoto 1987). 

9.8. Multi-Zone Models 

There are two multi-zone chemical evolution models for ellipticals available in the lit- 
erature (Martinelli et al. 1998; Tantalo et al. 1998). Martinelh et al. (1998) assumed 
that the elliptical galaxy is divided in several concentric shells of thickness Ai?^. The 
binding energy of the gas in each shell is computed after assuming a dark matter halo as 
described before. The model predicts that a galactic wind develops first in the external 
regions and then gradually in the more internal ones in agreement with the observed 
Mg2 vs. Vesc relation. As a consequence of this, the star formation lasts for a shorter 
time in the external regions with the consequence that an abundance gradient is created 
in agreement with observations. This model also predicts that we should observe higher 
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Figure 8. Predicted and observed color gradients for some field elliptical galaxies. From Menan- 
teau et al. (200f). Figure shows from left to right: 7gi4-band surface brightness map, V606-7si4 
color pixel map and V606-^8i4('") color gradient. Open circles represent observed gradients while 
solid lines are the model predictions obtained by means of the Martinelli et al. (1998) model. 



[Mg/Fe] ratios at larger galactocentric distances, the contrary of what should happens 
in the disks of spirals. In this case, we can speak of outside-in formation. Unfortunately, 
the available data inside galaxies do not allow yet to observe such an effect. It is worth 
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noting that this model can reproduce very well the observed color gradients in ellipticals, 
as shown in Figure 8. 

10. Conclusions on Ellipticals 

The main conclusions that we can draw from comparing theoretical models and obser- 
vations for ellipticals can be summarized as follows: 

• In order to explain the observed < [Mg/Fe] >*> in giant ellipticals the dominant 
stellar population should have formed on a time scale no longer than 3-5 -10^ yr, which 
corresponds to the time at which the SNe la rate reaches a maximum in these systems 
with strong star formation. 

• Uncertainties in the stellar yields of Fe and Mg from different authors can change 
the value of the [Mg/Fe] ratio but do not affect the conclusion above. 

• This observational finding argues against a hierarchical clustering formation scenario 
and favors the fact that ellipticals, especially those in clusters, are mostly old systems 
(see also Menanteau et al. 2001). 

• The increase of the [Mg/Fe] ratio with galactic mass suggests either that more mas- 
sive ellipticals are older systems or that the IMF is not constant among ellipticals or 
both. 

• Abundance gradients in ellipticals can be produced by biased winds and this would 
imply that [Mg/Fc] increases with increasing galacto centric distance. 

• Better calibrations for metallicity indices are necessary, especially taking into ac- 
count non-solar ratios in stellar tracks, before drawing firm conclusions. 

11. Evolution of Dwarf Galaxies 

Dwarf Irregular (DIG) and Blue Compact (BCG) galaxies are very interesting objects 
for studying galaxy evolution since they are relatively unevolved objects (see Kunth & 
Ostlin, 2000, for a recent exhaustive review). In bottom- up cosmological scenarios they 
should be the first self- gravitating systems to form, thus they could also be important 
contributors to the population of systems giving rise to QSO-absorption lines at high 
redshift (DLAs). In general, they are rather simple objects with low metallicity and 
large gas content, suggesting that they are either young or have undergone discontinuous 
star formation activity (bursts). An important characteristic of these systems is that 
they show a distinctive spread in their physical properties, such as chemical abundances 
versus fraction of gas. Matteucci and Chiosi (1983) were among the first in studying the 
chemical evolution of dwarf galaxies. They adopted analytical models as those described 
in section 3 and showed that closed-box models cannot account for the Z-log/U distribution 
even if the number of bursts varies from galaxy to galaxy, and suggested possible solutions 
to explain the observed spread: 

• a. different IMFs 

• b. different amounts of galactic wind 

• c. different amounts of infall 

Later on, Matteucci and Tosi (1985) presented a numerical model where galactic winds 
powered by SNe were taken into account. They concluded that different wind rates from 
galaxy to galaxy could explain the observed spread in 0,N vs. log/i but not the spread in 
the N/0 vs. 0/H diagram, suggesting that additional processes could have contributed 
to that. For example, different amounts of primary N from galaxy to galaxy. Kumai 
and Tosa (1992) suggested that different fractions of dark matter in different objects 
could explain the observed spread in the Z-log/U diagram. Pilyugin (1993) forwarded 
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the idea that the spread in the properties of these galaxies (i.e. He/H vs. 0/H and 
N/0 vs. 0/H) are due to self-pollution of the HII regions coupled with "enriched" or 
"differential" galactic winds. Differential winds should carry out of the galaxy certain 
elements more than others, for example the metals ejected by SN II should be favored. In 
more recent models (Marconi et al. 1994; Bradamante et al. 1998), the novelty was the 
contribution to the chemical enrichment of SNe of different type (II, la and lb) together 
with differential winds. In these papers the assumption was made that the products of 
SN II are lost more easily than the products of stars ending as WDs, such as C and 
N. Larsen et al. (2001) studied the chemical evolution of gas rich dwarf galaxies and 
concluded that primary N production from massive stars is not necessary to reproduce 
the N/0 vs 0/H and that the spread in this relation is likely to be due to the time-delay 
effect in the production of N relative to the production of O. In fact, in a starbursting 
regime when the starburst fades the elements produced by Type II SNe are no more 
produced whereas the elements produced on long timescales by single stars and Type 
la SNe are still produced. This causes the typical saw-tooth behaviour (see Figure 9). 
Therefore, the spread can be due to the fact that some objects are observed at different 
stages of the burst /inter hurts regime. They also concluded that ordinary winds are better 
than enriched ones in reproducing the properties of these objects. The existence of a 
luminosity-metallicity relation (although with spread) can be an indication for galactic 
winds acting more efficiently in low mass than in high mass high potential well objects. 

11.1. Evidences for Galactic Winds 

Meurcr et al (1992), Papaderos et al. (1994), Lequeux et al. (1995) and Marlowe et al. 
(1995) all suggested the existence of galactic winds in dwarf starbursting galaxies. The 
evidence is gathered by the indication of outflowing material travelling at a speed larger 
than the assumed mass of the objects. Papaderos et al. estimated a galactic wind flowing 
at a velocity of 1320 Km/sec for VIIZw403. The escape velocity estimated for this galaxy 
being 50 Km/sec. Lequeux ct al.(1995) suggested a galactic wind in Haro2=MKn33 
flowing at a velocity of ~ 200Km/sec. Martin (1996;1998) found supershells in 12 
dwarfs including IZwlS and concluded that they imply an outflow which in some cases 
can become a wind (namely the material is lost from the galaxy). 

Recent chemo-dynamical simulations for one instantaneous starburst also suggest the 
possibiliy of galactic winds and that these winds are metal enriched (MacLow & Ferrara 
1999; Recchi et al. 2001). 

11.2. Results for BCG from chemical models 

Purely chemical models (no dynamics) have been computed by several authors by varying 
the number of bursts, the time of occurrence of bursts, tburst, the star formation efficiency, 
the type of galactic wind, the IMF and the nucleosynthesis prescriptions (Marconi et al. 
1994; Kunth et al.l995; Bradamante et al.l998). The main conclusions of these papers 
can be summarized as follows: 

• The number of bursts should be N^ursts < 10, the star formation efficiency should 
vary from 0.1 to 0.7 Gyr''^ for either Salpeter or Scalo (1986) IMF but Salpeter IMF is 
favored. 

• Enriched winds, carrying out material at a rate proportional to the star formation 
rate, seem to be preferred (but see Larsen et al. 2001). The observed scatter in the 
observational properties can be due either to the winds or to the delay in the production 
of different elements. 

• If the burst duration is relatively short (no more than 100 Myr), SNe II dominate 
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Figure 9. The log(N/0) vs. 12 +log(0/H) for a sample of BCG. The data are from Recchi 
(2002). Overimposed are three models with a single burst of star formation and different star 
formation efficiency. In particular, the dotted line corresponds to an efficiency v = \Gyr~^ , the 
continuous line to i/ = 2.5Gyr~^ and the dashed line to v = 5Gyr~^. The burst duration is 100 
Myr. As one can see, the saw-tooth behaviour typical of a bursting mode of star fotmation is 
evident. 

the chemical evolution and energetics of starburst galaxies, while stellar winds seem to 
be negligible after the onset of SNe IL 

• The [0/Fe] ratios tend to be overabundant due to the predominance of Type II SNe 
during the bursts. Models with a large number of bursts Af6„rst=10 - 15 can give negative 
[0/Fe]. 

11.3. Results from chemo- dynamical models 

Recent chemo-dynamical models (Recchi et al. 2001;2002) assuming an instantaneous 
starburst but following in great detail the evolution of several chemical elements (H, He, 
C, N, O, Mg, Si and Fe) and adopting an efficiency for SN II energy transfer 775^/7=0. 003 
and for SN la /ysjv/o=l-0, suggest that the starburst triggers indeed a galactic wind (see 
Figure 10). In particular, the metals leave the galaxy more easily than the unprocessed 
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Figure 10. Chemo-dynamical simulation of an instantaneous starburst from Recchi et al. 
(2001). The development of a galactic wind is evident along the z-axis, since the assumed 
configuration is flattened. The age of the burst is indicated in Myr. 

gas and among the metals the SN la ejecta leave the galaxy more easily than the SN 
II ejecta. This is due to the assumed efficiencies for energy transfer for the two types 
of SNe. This assumption is in turn based on the fact that Type la SNe explode in an 
already heated and rarified medium (thanks to the Type II SNe which explode first) and 
therefore can transfer all of their energy into the ISM. As a consequence of this type of 
evolution the following conclusions can be drawn: 

• A selective loss of metals seems to occur in dwarf gas-rich galaxies. 

• As a consequence of the selective winds, the [a/Fe] ratios inside the galaxy are 
predicted to be larger than the [a/Fe] ratios outside the galaxy. In fact, the products of 
SNe la are lost more efficiently than those of SN II. 

• At variance with previous studies, most of the metals are already in the cold gas 
phase after 8-10 Myr owing to the fact that the superbubble does not break immediately 
(the SNe II inject only a fraction of their initial blast wave energy into the ISM) and 
thermal conduction can act efficiently. 
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• The model well reproduces the properties of IZwlS (the most metal poor galaxy 
known locally) if two bursts are assumed and are separated by 300 Myr interval. 

11.4. Dwarf galaxies and DLA Systems 

Finally, before concluding this section, we would like to draw the attention upon the 
fact that there are similarities between BCG, DIG and DLAs or more in general between 
DLAs and systems with a low level of star formation. 

The nature of DLA systems is under debate and the abundance ratios measured there 
can be used as a diagnostic to infer their nature and age. Matteucci et al. (1997) 
suggested that some DLA showing high N/0 ratio could be BCG suffering selective 
winds. In this respect, the similarity between DLAs and IZwlS may suggest that IZwlS 
is a survivor proto-galaxy which has just started forming stars. If this is true, dwarf 
irregular galaxies should be born at any time during the age of the universe, unlike 
elliptical galaxies which appear to have formed a long time ago and in a short time 
interval (see section 9). 

Plots of [a/Fe] vs. [Fe/H] and plots of [a/Fe] vs. redshift should be used to infer the 
nature and the age of these objects, when compared with chemical evolution predictions. 

The observed abundance ratios seem to indicate that DLAs show almost solar [a/Fe] 
ratios at low [Fe/H] (Pettini et al. 1999; Centurion et al. 2000). This is an indication 
that they are objects where the star formation proceeded slowly and that they have 
probably started to form stars long before the redshift at which we observe them. This 
could indicate that we are looking at the external regions of disks or at dwarf starbursting 
systems in the interburst phases (see Figure 11). However, more data are necessary to 
assess this point. One common problem related to the measurement of the abundances 
in DLAs is that some elements are dust depleted such as Si and Fe and therefore one 
should try to observe non-refractory elements such as N, O, S and Zn. One should also 
remember that Si and Ca do not strongly behave as a-elements since they are produced 
in a non-negligible way in SNe la. 

12. Chemical Enrichment of the ICM 

After having discussed the chemical evolution of galaxies it is important to conclude 
by studying how the evolution of galaxies can affect the chemical enrichment of the 
intracluster (ICM) and intergalactic medium. In the past years a great deal of work has 
been presented on the subject. The first work on chemical enrichment of the ICM was by 
Gunn & Gott (1972), Larson & Dinerstein (1975), Vigroux (1977), Himmcs & Biermann 
(1988). In the following years, Matteucci & Vettolani (1988) started a more detailed 
approach to the problem followed by David et al. (1991), Arnaud et al. (1992), Renzini 
et al. (1993), and many others. The majority of these papers assumed that galactic winds 
(mainly from ellipticals) are responsible for the ICM chemical enrichment. Alternatively, 
the abundances in the ICM could be due to ram pressure stripping (Himmes & Biermann 
1988) or to pre-galactic Population HI stars (White & Rees 197S). 

12.1. Models for the ICM 

Here we will describe briefly the metodology ckn'eloped by Matteucci & Vettolani (1988) 
(hereafter MV88). Starting from SN driven galactic wind models for ellipticals (see 
section 9) they computed the ejected masses vs. final baryonic total galactic masses for 
galaxies of different initial luminous masses (from 10^ to lO^^M©): 



(12.67) 




Figure 11. Rcdshift evolution of [Zn/H] for high surface brightness spirals (HSB) (bottom 
panels) and low surface brightness spirals (LSB) (upper panels) for different galactocentric 
distances in the disk, a) The evolution in metallicity of LSB formed at z=4 fits well the DLA 
values measured by Pettini et al. (1999).b) If LSB disks formed late, however, then they do 
not fit the data, c) Conversely, HSB disks that formed at z=4 become metal rich too quickly 
to explain the observations; d) only if there is a continuously forming population of HSB disks 
between z=4 and z=l then they can account for the metallicities of the DLAs. 



where i refers to either a single chemical element or the total gas. Then, they integrated 
M^-' over the cluster mass function obtained from the Schechter (1986) luminosity func- 
tion (LF) by assuming that only E and SO galaxies contribute to the chemical enrichment. 
This assumption was later confirmed by data of Arnaud et al. (1992). The total masses 
ejected into the ICM in the form of single species i and total gas are: 

Mtiius{> Mf) = Eifn*ih^kfnO-°-''^*^^'-^-^^^-T[{a+l+f3i), (M;/iV^)10"°-^*'^^"^-**^] 

(12.68) 

where a is the slope of the LF, / is fraction of ellipticals plus SO, M* is the mass at the 
"break" of the LF and is the magnitude at the break. 
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Table 1. Predicted and observed quantities 



Cluster 


(MFe)pred 


{MFe)obs 




{Mgas)ohs 


Perseus 


2.610" 


1.910" 


2.010" 


3.010" 


A2199 


1.310" 


1.210" 


10" 


1.510" 


Coma 


2.310" 


3.110" 


1.810" 


1.510" 


Virgo 


1.410^° 


1.610^° 


10^2 


2.010" 



Table 2. Predicted Fe abundances and abundance ratios 



Cluster Xfs/Xf,^ [Mg/Fe] [Si/Fe] 



Perseus 0.65 -0.65 -0.080 

A2199 0.65 -0.60 -0.096 

Coma 0.43 -0.63 -0.110 

Virgo 0.53 -0.64 -0.090 



12.2. MV88 Results 

MV88 considered 4 galaxy clusters: Perseus, A2199, Coma and Virgo for which /, n*, 
M* , Mpe and Mgas were known. In Table 1 we show the results they obtained for 
the 4 clusters concerning the total masses ejected by all galaxies in the clusters in the 
form of total gas and Fc, compared with the observed ones. It is immediate to see from 
Table 1 that their model could reproduce well the total amount of Fe but they failed 
in reproducing the total gas mass. However, this is not a failure of the model but the 
indication that most of the gas in clusters has a primordial origin, namely has never been 
processed inside stars. The same conclusion was reached later by David et al. (1991) 
and Renzini et al. (1993) among others, but see Chiosi (2000). 

Therefore, the Xpe/ Xp^^,^ shown in Table 2 was calculated as {Mpe)pred/ {Mgas)obs- 
The predicted Fe abundance in the ICM relative to the Sun is in agreement with the 
observations then and now {Xpg/{XFe^)obs = 0.3 — 0.5 (Rothenflug & Arnaud 1985; 
White 2000). Low values for [Mg/Fe] and [Si/Fe] were predicted due to the assumption 
that all the Fe, produced by Type la SNe, was soon or later ejected into the ICM, as 
shown in Table 2. With Salpeter IMF, Type la SNe contribute > 50% of the total Fe. 

12.3. [a/Fe] Ratios m the ICM 

Models including Type la and II SNe predict an asymmetry in the [q/Fc] ratios (> 
inside the ellipticals and < in the ICM) due to the different roles of SNe II and la in 
Fe production (Renzini et al. 1993.) ASCA results (Mutshotzky et al. 1996) originally 
suggested [a/Ye]jcM > (-1-0.2 dex). However, Ishimaru & Arimoto (1999) pointed out 
that [a /Fc] ~ in the ICM if the metcoritic Fe abundance is adopted instead of the 
photospheric value adopted in the other paper. On the basis of the ASCA results on the 
overabundance of the a-elements, Matteucci & Gibson (1995) discussed how to reproduce 
the [a/Fe]>0 ratios both in stars and ICM. They adopted the same SN feedback as MV88 
and dark matter halos were also included in the galaxy models. They concluded that 
it is possible to obtain overabundances of the a-elements in the ICM only if not all of 
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the produced Fe is lost from the galaxies (i.e. only early winds). On the other hand, 
MV88 had assumed that all the Fe is soon or later ejected into the ICM. They concluded 
that a flat IMF (x=0.95) plus only early winds can reproduce [a/Fe]icM > 0. A similar 
conclusion was reached by Loewenstein & Mutshotzky (1996). However, the situation is 
not very realistic since the gas in cluster galaxies is likely to be stripped soon or later 
because of environmental effects, and the [a/Fc] ratios in the ICM are likely to be solar or 
undersolar. More recently, Martinelli et al. (2000) recomputed the ICM enrichment by 
adopting a more realistic model for the evolution of ellipticals (multi-zone). The model 
(already described in section 9) predicts a more extended period of galactic wind and 
more metals and gas in the ICM than the one-zone model with only early winds. They 
computed the evolution of abundances vs. redshift for a constant LF and concluded 
that there is no evolution between z=l and z=0 and [a/Fe]/cM <0 at the present time. 
Very recently Pipino et al. (2002) computed the chemical enrichment of the ICM as a 
function of redshift by considering the evolution of the cluster luminosity function and 
confirmed the conclusions of Martinelli et al. (2000). Very recent data from XMM- 
Newton (Gastaldello & Molendi, 2002) seem to indicate a low [0/Fe]< ratio for the 
ICM, thus supporting models where Type la SNe are efficient and all the produced Fe is 
ejected into the ICM. In Figure 12 we show the recent results of Pipino et al. (2002). In 
the upper panel of Fig. 12 is indicated the predicted evolution as a function of redshift 
of the thermal energy per particle in the ICM due to galactic winds from E and SO 
galaxies.In the lower panel is shown the evolution of the total mass of Fe ejected by the 
cluster galaxies into the ICM as a function of redshift. It is worth noting that 1 keV per 
particle is the energy required to explain the observed Lx — T relation in clusters (e.g. 
Borgani et al. 2001). The model in Figure 12 shows that the energy injected by SNe is 
not enough (at maximum 0.4 keV) and that other sources of energy are required, such as 
active galactic nuclei and QSO. Figure 12 also shows how important is the contribution 
from Type la SNe to the chemical enrichment and energetic content of the ICM. In fact, 
SNe II are important in triggering the galactic wind but after star formation stops the 
wind is sustained only by Type la SNe. 

12.4. [a/FeJ ratios and IMLR 

Abundance ratios and the Iron Mass to Light Ratio (IMLR) are good tests for the 
evolution of galaxies in clusters since they do not depend on the total cluster gas mass. 
The IMLR is defined as (Renzini et al. 1993): 

IMLR = M^p^^/Lb (12.69) 

and the IMLR oc For Hq = 50 the IMLR=0.02 (Mq/Lq) (Arnaud et al. 1992) 

practically constant among rich clusters but it drops for poor clusters and groups. On 
the one hand, the IMLR can impose constraints on the IMF in cluster galaxies, on the 
baryonic history of the clusters/groups and on the number of SNe exploded in galaxies. 
In fact, the most straightforward interpretation of the constancy of the IMLR among 
rich clusters is that they did not loose Fe at variance with the groups or small clusters 
having a lower IMLR indicating a loss of baryons (Renzini, 1997). On the other hand, 
the [a/Fe] ratios can impose constraints on stellar nucleosynthesis, IMF, different roles 
of SNe and SN feedback. 

13. Conclusions on the ICM 

The main conclusions on the chemical enrichment of the ICM can be summarized as 
follows: 



36 Francesca Matteucci: Chemical Evolution of Galaxies and Intracluster Medium 



E elected in the ICM. 

pp 

produced by SNIa. 

E produced by SNII. 



redshift 



Mpj. (Mg) ejected in the ICM. 
Mp^ (Mq) produced by SNIa. 
Mp^ (Mq) produced by SNII. 



redshift 



Figure 12. The predicted evolution ol the thermal content of the ICM (upper panel) as a 
function of redshift. The different contributions from SNe of different type is indicated. In the 
lower panel is shown the evolution of the total Fe mass as a function of redshift. The assumed 
cosmological model is Q.m = 0.3, ^\ = 0.7 and h=0.70. 



• Good models for the chemical enrichment of the ICM should reproduce at the same 
time the [a/Fe] ratios inside galaxies and in the ICM. They should also reproduce the 
IMLR. 

• Crucial parameters are: SN feedback, SN nucleosynthesis, IMF and whether all the 
stellar ejecta (produced over a Hubble time) can reach the ICM or remain bound to the 
parent galaxy. 

• Type la SNe play a fundamental role in the chemical and energy content evolution 
of the ICM. 

• Good models for the ICM enrichment, assuming that the energy transfer from SN 
II is ^ 3% and from SNIa is 100%, predict an energy per particle in the ICM of Ep ^ 0.4 
keV, not enough to break the self- similarity in clusters. 
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